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We calculate the rate for pair production of light pseudoscalars by strong inhomogeneous and 
static electric (magnetic) fields. We show that, in the case of axions, the stability of atoms over 
cosmic lifetimes is jeopardised unless the Peccei-Quinn symmetry breaking scale /„ is larger than 
O(10 w GeV). 



The Schwinger mechanism is a non-perturbative pro- 
cess by which an infinite number of zero frequency pho- 
tons can produce an electron-positron pair. For exam- 
ple a constant electric field larger than 0(lQ ll V/cm) 
will decay into electron-positron pairs, according to the 
Schwinger pair production probability formula Q . In this 
paper we consider the pseudoscalar-two-photon coupling. 
We find that the pseudoscalar pair production probabil- 
ity is given by an expression w ~ g 2 exp —(0(m 2 /g)) in 
terms of the pseudoscalar-two-photon coupling g and the 
pseudoscalar mass m. This allows us to constrain the 
parameter space of {m, g} for any pseudoscalar model, 
based on the stability of electromagnetic fields. In ax- 
ion models [00] the effective axion-two-photon coupling 
arise from axion couplings to charged fermions in a one 
loop diagram. This coupling is inversely proportional to 
the Peccei-Quinn symmetry beaking scale f a , and is in- 
dependent of the mass of the fermion in the loop. The 
axion mass is inversely proportional to the square of f a . 
The pair production formula therefore can be expressed 
in terms of the single parameter f a . 

The generic electromagnetic interaction of a pseu- 
doscalar (f> can be written as (see Fig.l) 



Cj = \g4>F^F^ 



(1) 



with F M „ = dfiAi, 



duAf, and W v = \e^ af3 F a p. 
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FIG 1: Axion-two-photon interaction. 

We need to evaluate the loop diagram of the type 
shown in Fig. 2 with infinite number of zero-frequency 



photon external legs. The imaginary part of this dia- 
gram gives the probability for the decay of the external 
electromagnetic field. 






FIG 2: The effective photon-photon interaction diagram 
on integrating out the axion and a photon field. 

We evaluate this diagram in two stages. First we find 
the effective vertex for the two axion-two photon vertex 
which arises when we contract the photon legs between 
two vertices like (fil). This two axion-two photon vertex 
is of the form 



\g<i>) M r/ '"- 



~Wv 



k 2 



■i-k^F" 



(2) 



The factor of 4 in equation (0) is for the four possible 
ways of joining the photon legs. Due to the presence of 
the fc 2 term in the denominator, (0) is non-local. How- 
ever, momentum k is integrated over when we calculate 
the effective action for the external electromagnetic field. 
One can therefore make use of the identity 



d A k f{kji 



IM^IJ,' 



d 4 kf 



\ 4 



to simplify (El) and obtain 



C'i = \g 2 <j?F, v F^ = i 5 V(E 2 - B 2 ; 



(3) 



(4) 



When the momentum k is integrated over, the effective 
two axion-two photon interaction reduces to a local in- 
teraction vertex (Fig 3) given by (||). 
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FIG 3: The effective two axion-two photon interaction 
in the loop. 

Next, to calculate the diagram with infinite insertions of 
the two-photon vertices (0) in an axion loop we use the 
Schwinger method. Starting with the total action for the 
axion-photon system we integrate out the axion fields to 
determine the effective action for just the electromagnetic 
fields 

e »S e/ /[E,B] = [ x>4> e*-/ d4x {i<M- 92 - m2 )*+29V( E2 - B2 )} 

To proceed we shall evaluate £ e // using a technique de- 
veloped in S. It amounts to expanding the classical 
background fields about a reference point, to compute 
the (interacting) Green's function in momentum space, 
to revert it to configuration space, and finally obtaining 
Ceff by a parametric integration. We shall be inter- 
ested in stationary, inhomogeneous (on the scale of the 
pseudoscalar Compton wavelength) external fields (elec- 
tric and/or magnetic fields). The Green's function of the 
axion field in a background electromagnetic field is given 
by 

[8 2 + to 2 - a - pi{x - x'Y - 7 2 (x - x')\x - x'y] 

xG(x,x') = 5 4 (x-x') (5) 



[i, j = 1,2,3) where, as advertised, we have expanded the 
electromagnetic fields to second order about the reference 
point x'. In momentum space this equation reads, 







d 2 



-Po + v\ + m 2 -a + ifa- h 7y ., 

dpi J dpidpj 



G(p) = l (6) 



This equation is solved by making the ansatz 

G{p) =i ds e ^(Po-™ 2 ) e -*P-A-P+B. P +c ^ 

and solving for matrix A = A y , vector B = B l , and C. 
The solution is given by 

A = — *y • tanh 27s 

B = -\l~ 2 ■ (1 - scch27s) • (3 

C = ias tr In cosh 27s 

2 ' 

+ ~ j3 ■ j- 3 • (tanh 2 7 s - 2 7 s) • (3 (8) 
8 

where tr is a trace over indices i,j. 
Since 



dC 



eff 



da 



— G(x,x) 

d A p 
2i J (2tt) 4 



G(p) 



(9) 



subject to the boundary condition C e ff = when E = 
B = 0, we can find C e ff by gaussian integration to obtain 



-eff 



32tt 2 



00 ds 

~3 



- lS (m 2 -a) / det 



27S 



sinh 27s 



1/2 



Jl(s) 



(10) 



where 



i{s) = ~p-r 



(tanh 7s — 7s) • (3 



We can now apply this formula to the special case of 
spherically symmetric electric fields, like the ones in 
atoms. In this case g 2 E 2 appearing in the interaction 
lagrangean is written as, 



g 2 E 2 



9 E Q 



2 j^2l f 

9 E (r 



r*) + \ 9 2 E 2 "{r 



ro 



(11) 



when expanded about ro. Clearly, a — g E 2 . Also, 

n 2 E 2 ' 

vector /3 l above reads now (3 = - — 2-ro and matrix 
7 2 is, after diagonalization, j 2 = diag(7 2 ,7 2 ,7 2 ) where 



= \g 2 E 2 " and 7 f = 



IK 

2rn 



tion l(s) given in equation(l2) is explicitly 

l(s) = rl^li 3 (tanh 71s - 71s) 
So, finally, 



With this input the func- 



(12) 



£eff - --■ 



i r°° ds 

32tt 2 J s 3 


c - is (m 2 - a ) 2 72S 

sinh272S \ 


1 2 7is c a(s) 
j sinh27is 


— ism" 



(13) 



Now, when we specialise to the field created by an 



atomic nucleus E 2 — IMh we have, 7? 



2 _ g Z z a 



27T,'!; 



and 



7i 



2 _ 5g 2 Z 2 a 



2irr° 



The probability w per unit volume and unit 
time for creation of pseudoscalar pairs is 2ImC e ff and 
can be obtained from the previous equation by complex 
integration. We obtain: 



^^(-i)"a ie - K ( 2 " +1 )- 

n=0 



Here, 72 = +(— 7 2 ) 2 7 an d k 

A = 72 4 7r 3 ^o- The 
series above are given by 



K^ A )- 



(14) 



with 



A = 72 7! 3 rg. The coefficients C n in the asymptotic 



t ..-.„ 



du ■ 



^—(ku+\ cot 77 ) 



[u+ (2n + l)ir} 2 (sinu) 2 sinh 22- 



+ (2n+l)7r] 



(15) 



For the case of axions, the mass and coupling are related 
to the global symmetry breaking scale f a [0J|] 



<?~0.75 7 



2irfa 

10 7 GeV 
m a ~ 0.6 eV 

J a 



(16) 



and one can estimate the magnitude of the pair-creation 
phenomenon in terms of the single parameter f a . The 
asymptotic expansion in the equation for w only makes 
sense for k > 1 and w is vanishingly small for k> 1, 
i.e. for second derivatives of the fields small on the scale 
of m a . Only when fields approach critical values for 
which k ~ 0(1), pah creation could eventually reach 
catastrophic rates and produce vacuum breakdown and 
disruption of the external field. Pair creation, which is 
driven by the second derivative of the field, becomes crit- 
ical (i.e. k = 1) in a Z-atom for r cr n — ( ^ O 1 ^ 3 frn. 
Since r crit must be larger than a few /m, our arguments 
apply only for 



f a >O(10 2 GeV). 



(17) 



Let us now constrain f a by realising that the observed 
cosmic abundance of atomic helium agrees with nucle- 
osynthesis up to a few percentage and therefore we can 
use the stability over the lifetime of the Universe of He- 
atoms. We consider a thin shell of thickness Sr = 0.1 r 
at critical radius r crit inside the He-atom. The atom will 
become ionised due to breakdown of the electric field in 
this shell at some time r when the product Vwt ~ 1 
where V is the volume of the shell, and w is as given in 



Imposing r > 1.5 x 10 10 years, and taking into account 
(J17|), we get the following exclusion range for /„: 



10 2 GeV < f a < 0.6 x 10 10 GeV 



(18) 



We note that, for the whole of the exclusion region, 
the critical radius r cr i t <S 1A, i.e. inside the atom, and 
the Compton wavelength of the axion exceeds 10 3 fm, 
much larger than the distances over which the field is 
inhomogeneous . 

We can improve the constraint on f a by recognising 
that there are many heavy atoms which have remained 
un-ionised over the lifetime of the universe. We know 
from astronomical observations that there are many stars 
with lifetime larger than r ~ 10 10 years. This is deter- 
mined by comparing the abundance of heavy elements 
ratio (Th/Eu) in the stars with that in the solar system. 
Using the fact that the lifetime of thorium is 14 x 10 9 
years it has been estimated that there are stars which 
have ages in excess of 12 x 10 9 years ||. Since these 
atomic abundances in the stars are measured from the 
intensity of their spectral lines we know that they have 
not been ionised due to axion-pair production over times 
larger than r ~ 10 10 years. However, because r cr i t grows 
with Z, we cannot use a too large a Z for then r cr n 
would penetrate the atomic electron cloud and screening 
would distort the field. We estimate that we can go up 
to Z ~ 20 and then the bound can be pushed up only 
to somewhat above 10 10 GeV . Anyway, we are not aim- 
ing at a very elaborate limit here M. We just wish to 
point out that to avoid the breakdown of electric fields 
in atoms by catastrophic axion pair-creation requires f a 



to lie about or beyond the 10 10 GeV scale, which is in 
the ballpark or already closing the presently still allowed 
axion-photon coupling window. 
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